3D Jet Tomography of the Twisted Color Glass Condensate 
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Jet Tomography is proposed as a new test of Color Glass Condensate (CGC) initial conditions 
in non-central A + A collisions. The kr factorized CGC formalism is used to calculate the rapidity 
twist in the reaction plane of both the bulk low pr < 2 GeV matter as well as the rare high pr > 6 
GeV partons. Unlike conventional perturbative QCD, the initial high pr CGC gluons are shown to 
be twisted even further away from the beam axis than the the low pr bulk at high rapidities \n\ > 2. 
Differential directed flow vi(pr > 6, \n\ > 2) is proposed to test this novel high pr rapidity twist 
predicted by the CGC model. 

PACS numbers: 12.38.Mh; 24.85.+p; 25.75.-q 



Introduction: It was pointed out in Ref.Q] that the 
QCD matter produced in high energy noncentral A + A 
nuclear reactions violates (locally) Bjorken longitudinal 
boost invariance in the transverse plane even if the global 
rapidity distribution, dN/drj, is independent of r\ = 
sintC (p z /m±). The intrinsic longitudinal boost non- 
invariance occurs even in symmetric A + A reations be- 
cause locally in the transverse xj_ plane, the initial gluon 
density, p g (?7,xj_) = dN g /dr]d 2 x.±, has a generic "trape- 
zoidal" form in the rapidity variable. This peculiar struc- 
ture arises in noncentral (b > 0) collisions because there 
is an asymmetry between the local number of interacting 
projectile and target nucleons, A^(x^;6 > 0) ~ A 1 / 3 , 
that can vary by an order of magnitude in heavy nu- 
clei with the transverse coordinate, xj_. Since low pr 
bulk matter is known to be produced proportional to the 
number of participating nucleons, the slope of the rapid- 
ity trapezoid varies with xj_ as dp/dr/ oc Ai/(x_l; b)/2Y ~ 
A 1 / 3 / logs, where 2Y = logs is the rapidity gap between 
the projectile and target nuclei. At infinite energies this 
slope approaches zero, but at the Relativistic Heavy Ion 
Collider (RHIC) and the Large Hadron Collider (LHC) 
it may be large enough to be observable via 3D (p± , 4>, rf) 
extensions^ of jet tomography^- 

The trapezoidal rapidity asymmetry has been ob- 
served at all energies in p + A reactions and it was also 
clearly seen in D + Au reactions at RHIC|3j, |4|. The 
data are well reproduced by phenomenological soft+hard 
(string+mini-jet) models such as 0,0,0- However, these 
trapezoidal features are also well reproduced by gluon 
saturation models such as in the Kharzeev-Levin-Nardi 
(KLN)H implementation of the Color Glass Condenstate 
(CGC) theory |3, llCj of low x parton initial conditions. 

However, the phenomenological soft+hard and CGC 
approaches differ significantly in their predictions for 
moderately high pt partons above Q s . In the former 



case, the high pt parton production is calculated from 
the collinear factorized on shell gg — > gg approxima- 
tion of pQCD for hard processes. In the KLN/CGC ap- 
proach both soft and hard partons are calculated using 
the kx factorized off-shell gg — ► q gluon fusion approxi- 
mation to QCD parton production[ll|. The advantage of 
collinear factorized approximations is that they employ 
experimentally well determined nucleon parton structure 
functions, xGn(x,Q 2 ). In addition, corrections beyond 
lowest order can be systematically evaluated. The main 
disadvantage in applications to A+A is that possible non- 
linear nuclear effects at low x could significantly modify 
the assumed linear relation G a = AGn in different kine- 
matic regimes. 

In contrast, KLN/CGC predictions are based on the 
convolution of unintegrated nuclear gluon distributions, 
(f>A{x,k^) via the Gribov-Levin-Ryskin (GLR) formula. 
This approach has the advantage of including a nonlin- 
ear gluon evolution as x — > via a single gluon satura- 
tion scale Q s (x, A) that can be determined in the weak 
coupling but strong field approximation. The main dis- 
advantage of this approach is the uncertainty related to 
how to take the A — * 1 nucleon limit, (f>j^(x,k^). This 
limit is important phenomenologically because nuclear 
modification is measured relative to the N + N baseline. 
In addition, near the nuclear surface the participant nu- 
cleon density decreases rapidly and Q s is driven below 
1 GeV where the weak coupling strong field approxima- 
tions become suspect. Since the surface regions in fi- 
nite A+A contribute significantly (~ 10 — 20 %) to the 
global dN/dy, the specific implementation of the low A 
limit is important to devise experimental tests sensitive 
to this limit of CGC. Unlike the integrated parton dis- 
tributions, there is no general consensus yet on the form 
of 4>n 12]. Additional theoretical uncertainty is asso- 
ciated with the unknown applicability range of the first 



order GLR kx factorization formulation [llj as com par ed 
to proposed higher order nonlinear generalizations 13]. 

At RHIC the strongest support for the KLN/CGC ap- 
proach is its remarkable ability to reproduce the exten- 
sive systematics of the energy and nuclear size depen- 
dence of the global px integrated dN^/dy. This results 
from a specific dependence of the saturation scale, Q s , 
on -y/s and A. In contrast, phenomenological soft-hard 
models such as HIJING fail to reproduce the systematics 
because the separation scale, po ~ 2 GeV, between soft 
and hard parton production was assumed to be indepen- 
dent of those variables. It is an important open question 
of how high in pt can the fcT-factorized KLN/CGC ap- 
proximation be pushed in specifying the A + A initial 
conditions versus how low in px can the conventional 
collinear factorized approximation be pushed. Both ex- 
perimental and theoretical control over the initial condi- 
tions in A + A at RHIC are essential to strengthen the 
current case for the discovery of new forms of matter, 
the strongly coupled Quark Gluon Plasma (sQGP) and 
CGC, at RHIC Hill [13. 

In this letter a new jet tomographic approach is pro- 
posed based on extending the discussion in Ref.[l| to 
KLN/CGC initial conditions of the sQGP bulk as well its 
extrapolation to high px Q s jet partons. The idea is to 
exploit the difference between the geometric distributions 
of jets relative to that of the bulk matter as illustrated in 
Fig^ We focus on the predicted azimuthal dependence 
of the jet quenching pattern, RAAiVi 4>iP±\ b > 0)j an d 
long range rapidity correlations induced by the generic 
intrinsic rapidity twist of the bulk matter. 

At midrapidity, r\ = 0, the elliptic asymmetry of the 
reaction geometry in noncentral (b > 0) in A + A re- 
actions leads to a well known elliptic asymmetry in jet 
quenching. However, the rapidity twist has no observable 
effect at mid-rapidity. At positive rapdities 77 > 2 , on 
the other hand the rapidity twist of the bulk shifts the 
center of mass away from x = 0, while at negative rapidi- 
ties that shift has opposite sign. In conventional pQCD, 
collinear factorized gg — > gg predicts a jet distribution 
that is proportional to the local binary collision density, 

cr NN T B (r+)T A (r-), where r± = \JJx± f) 2 +y 2 , and 
Ta(x) is the Glauber nuclear profile functional . The 
jets are therefore produced symmetrically about xj_ = 
at all 77. The collinear factorized jet density therefore has 
no rapidity twist as illustrated by the grey ellipse in Fig^ 
The bulk matter density is effectively rotated away 
from the beam axis because the bulk density varies ap- 
proximately as {{Y - r))T A (r + ) + (Y + r])T A (r^)}/2Y, 
which is not reflection symmetric about xj_ = away 
from midrapidity. At px < Q s the CGC model pro- 
duces approximately the same rapidity twist as wounded 
nucleon string models since this is a direct consequence 
of local participant versus binary collision scaling of the 
bulk. However, due to the nonlinear equation determin- 




FIG. 1: (Color Online) Illustration of the initially twisted 
sQGP gluon densityfil relative to the beam axis in the (x,rj) 
reaction plane. CGc||, BGK0, and HIJING predict this 
generic low px < Q s locally boost non-invariant structure. 
Also shown are the relative rotations of the high pr 2> Q B jet 
partons in the kx factorized CGC model as well as conven- 
tional collinear factored pQCD. Jet quenching through the 
sQGP leads to opposite sign first azimuthal harmonic mo- 
ment, vi(px S> Qs,rj) in the two approaches. The projectile 
and spectator nuclei are indicated by half circles together with 
the sign convention of low px directed flow v\ . 



ing the local saturation scale Q s (xj_,x) the bulk density 
surface region can sharpen significantly over conventional 
Wood-Saxon geometry included in Ta(xj_). This sur- 
prising change of the bulk surface geometry in the KLN 
implementation of CGC was first pointed out by Hirano 
and Naraflij. 

However, as we show below, at higher p T the CGC pre- 
dicts even greater rapdity twist away from the beam axis 
than the bulk as shown below and illustrated in Fig^ In 
this paper, we explore some tomographic consequences of 
using the kx factorized formalism and CGC type unin- 
tegrated gluon distributions to produce both the bulk 
and jet matter in the reaction. This anomalous rapid- 
ity twist effect is opposite to that discussed in This 
effect occurs because the different nuclei are probed at 
asymmetric Bjorken momentum fractions while produc- 
ing high px matter. 

The Local Gluon Distribution: The local generalization 
of fcT-factorization GLR formula [Tl| used by KLN0 and 
Hirano and Nara 0] is given by 

dN g = 4tt q.(p» ) ^ x 
dpxd 2 xxdrj Cf Pt J 

,kx+px 2 _ k T -pr n2 - \ m 

Up = 2 c Nc and the collinear momentum fractions are 
given by kinematics, X1.2 — px exp(±rf)/^/s. The QCD 
coupling, a s , is evaluated at p\ and regulated at low px 
by imposing a maximum value a max = 0.5. Note that 
we use 77 to denote the rapidity rather than the pseudo 
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rapidity as in Q. 

4>a,b are the unintegrated gluon distributions which, in 
principle, possess a Bjorken x dependence determine d by 
nonlinear evolution equations of the CGC theoryjjj ^| 
and their kx dependence is fixed by a characteristic satu- 
ration momentum, Q s (x) . In the McLerran-Venugopalan 
approach [l9| the gluon distribution is suppressed below 
the saturation scale 4>a ~ \og{Q 2 s /k T ) compared to the 



pcrturbative form 



We use a parameteriza- 



tion similar to the KLN model approach as used in [l8j |. 
However, we use (for numerical convenience) the follow- 
ing Lorentzian form of 4>a.b for all values of kx- 



4> A {x,k T \x T ) 



Q 2 



s,A 



Q s 



A 2 



(2) 



The momentum scale A = 0.2 GeV is a regulator for the 
high rapidity y > 4.5 region as in 18]. The constant 
k ~ 0.5 is a parameter set to reproduce dN g /drj ~ 1000 
at midrapidity central collisions. The transverse coordi- 
nate dependence is implicit in the saturation momentum 
determined numerically for each nucleus. 



v {x,x T ) 



2tt2 
C f 



as(Q 2 S! A) xG nuc(x,Q 2 SiA )T A (x T ), (3) 



where Ta is the Glauber profile of nucleus A. We use 
standard diffuse Woods-Saxon profiles. The projectile 
and target nucleii are set up such that the spectator v% 
is positive at forward rapidity. 

The KLN parametrization is used for the nucleonic 
gluon distribution. 

xG mc {x, Q 2 ) =Klog( Q + A )x~ x (l - x) n (4) 

A QCD 

The momentum scales A and Aqcd are set to 0.2 GeV. 
The x~ x term accounts for the rapid growth of small 
x gluons while the factor of (1 — x) n was introduced in 
KLN to account qualitatively for the rapid depletion of 
gluons as x — > 1 outside the small x framework of the 
CGC model. As in the KLN approach, we set A = 0.2 
and n = 4. K ~ 1.35 is used to set (Q 2 s (x = 0.01)) ~ 2 
GeV 2 for central collisions at midrapidity. 

Once the distribution shown in Eq. \I\is evaluated, we 
are ready to investigate the transverse plane dependence 
(and hence the intrinsic spatial twist) of the produced 
gluons as a function of px and r/. We measure this shift of 
material away from the centre of the transverse reaction 
plane by calculating the average horizontal transverse co- 
ordinate. 



(x)(j>t,v) 



J d XTxdNg/dpxd x T drj 
J d 2 xxdN g / ' dpxd?xxdrj 



(5) 



Hereafter, x is the transverse coordinate in the direction 
of the reaction plane, b. Positive x points toward the 
projectile rj = Y spectator displacement . Fig. [21 shows 
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FIG. 2: (Color Online) The figure shows the average trans- 
verse spatial coordinate (x) for produced gluonic matter in 
the Brodksy-Gunion-Kuhn [f| type twist model as in as 
well as in the current CGC model. Notice that, for BGK, 
there is no dependence on px other than the assumed lack of 
twist for matter with px > 3 GeV. The CGC model, however, 
has a pronounced greater twist for higher px > 6 GeV. 



(x) as a function of px and r\ for the hybrid Brodsky- 
Gunion-Kuhn (BGK) |5j participant and binary jet pro- 
duction model (as used in Q) as well as for the current 
CGC model (for b = 9 fm). The rapidity twist of the 
BGK model is seen by the increasing (x) as rj increases 
for the bulk px < 3 GeV matter. 

While the bulk (x) of the CGC model is similar to the 
bulk BGK shifts, the high px shifts behave oppositely. 
In the two component BGK approach (x) = for high 
Pt- We show in Fig. [3 that the rapidity twist (d(x)/dr]) 
increases at high px > 6 GeV. 

Tomography and the Inverse Twist: Jet tomographic 
analysis uses the attenuation of jet matter while passing 
through the bulk in order to gain information about the 
density profile of the bulk . The observable used most 
commonly in tomographic analysis is the nuclear modifi- 
cation factor, Raa, which measures the deviation of the 
produced nucleus-nucleus spectrum, if any, from a simple 
binary scaled p-p spectrum. The twist effect investigated 
in the previous section can be observed by looking at the 
Raa{pt, Vi 4*) °f jets in the transverse plane. 

The azimuthal dependence of Raa will change as a 
function of rj for a given px jet due to the differing twist 
of the jet distribution over 77. In order to calculate the 
Raa{pt ■JJjA') we use the geometric model of Drees, Feng 
and Jia |2fJ. The nuclear modification factor is obtained 
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FIG. 3: (Color Online) The directed flow v\ as a function of 77 
for different p T at b = 9 fm. Both the CGC model and BGK 
model are given for comparison. 



by, 



Raa(pt,ii, 



J cPxre 



-fix(zT,4>,v) ^N g 

dpTd' 2 XTdr] 



Jd 2 X T 



dp-j<d' 2 XTdri 



(Pt,v) 



(6) 

/i = 0.04 is the parameter used to set Raa(v = 0, b = 
0) ~ 0.25. Opacity, x, is the line integral over the bulk 
distribution that is the cause of the attenuation experi- 
enced by the jet, calculated as in j^j. The length de- 
pendence of opacity is characteristic of radiative parton 
energy loss in Bjorken expanding matter. 

We can use Eq. [H] to probe the twist. The anti-twist 
effect can most easily observed via the first azimuthal 
fourier moment of Raa(pt,V> 0), the directed flow v%. 
Fig. |2| shows vi as a function of r\ for different values of 
the pt- Note that for all px values there exists a rapidity 
at which the directed flow flips sign. This flip occurs 
at lower values of the rapidity for higher values of the 
Pt- The sign of v\ is positive in CGC at high rapidity 
in the same direction as the projective spectator. The 
change in sign is a novel prediction using the KLN/CGC 
model. In conventional factorized QCD jet production, 
the high px v\ is negative and in the same direction as 
the low pt bulk directed flow but increasing with px as in 
hydrodynamics. Both models exhibit long range rapidity 
anticorrelations of the v\ . This anticorrelation may make 
it easier to measure the sign of v\ and test up to how high 
Pt can the KLN/CGC model describe initial conditions 
at RHIC. 
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